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Abstract 

The constitutive quantities in Mori's theory, the residual forces, are expanded in 
terms of time dependent correlation functions and products of operators at t = 0, 
where it is assumed that the time derivatives of the observables are given by prod- 
ucts of them. As a first consequence the Heisenberg dynamics of the observables 
are obtained as an expansion of the same type. The dynamic equations for corre- 
lation functions result to be selfconsistent nonlinear equations of the type known 
from mode-mode coupling approximations. The approach yields a neccessary con- 
dition for the validity of the presented equations. As a third consequence the static 
correlations can be calculated from fluctuation-dissipation theorems, if the observ- 
ables obey a Lie algebra. For a simple spin model the convergence of the expansion 
is studied. As a further test, dynamic and static correlations are calculated for a 
Heisenberg ferromagnet at low temperatures, where the results are compared to 
those of a Holstein Primakoff treatment. 
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1 Introduction 



In Mori's theory [1] the dynamics of a set of relevant observables in the Heisenberg 
picture are transformed to an equation of motion of the Langevin type showing a 
systematic part and a residual (stochastic) force as it is suggested by the phenomeno- 
logical theory of Brownian motion. The great success of Mori's theory is due to the 
fact that the rewritten equations of motion allow for excellent approximations in 
calculating linear response functions which are determined by the systematic part of 
Mori's equations. 

In treating the systematic part of the equation the main problem is to get an adequate 
approximation for the integral kernel given by the correlation functions of the residual 
forces. For its evaluation mainly two methods have been very successful: A simple 
perturbation theoretical calculation (for an introduction see [2]) and a factorization 
procedure into correlation functions of the observables (mode-mode coupling, e.g. 
[3,4] and references therein). 

From a general point of view, the basic quantity in Mori's theory is the residual 
force. The purpose of the present paper is to point out that for a large class of 
systems of interacting particles or spins, one can find a systematic approach for the 
residual force itself, thus being able to go beyond the calculation of linear response 
functions in terms of static correlation functions: From the residual force one can 
deduce approximations for 

i) the Heisenberg dynamics of the relevant operators 

ii) the dynamic correlation functions of the observables 

iii) the static correlations 

We will show that for systems as interacting spins, where the time derivatives of the 
relevant observables G are given by a superposition of products of the observables 

[H, G\ = V {1) G + V {2) GG + V^GGG + ■■■ (1.1) 

one can expand the residual force in terms of the dynamic correlation functions. This 
leads to an approximation scheme for the points (i) and (ii), and together with a 
Lie algebra for the observables to point (iii). In using this scheme one can see that 
the calculated Heisenberg dynamics do not have divergent secular terms which occur 
in a simple perturbation treatment. The dynamics of the correlation functions are 
governed by a set of nonlinear equations which are close to the usual mode-mode 
coupling equations. The static correlations in these equations can be expressed by 
dissipation-fluctuation relations and the Heisenberg dynamics which are given by 
the dynamic correlation functions and operators at time t — 0. Thus all quantities in 
principle can be calculated from the bare interactions V^ 2 \ 
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As already mentioned the essential condition for our treatment is that the time 
derivatives [H, G] of the relevant observables are given by the products of the ob- 
servables (Eq.(l.l)). Additionally we exploit the unitarity of the time evolution in 
Hilbert space which written with the Liouvillian 



is not automatically included in a Liouville space calculation, since G and GG are 
just two different elements of the space. 

It is to be expected that our expansion scheme could be generalized to the case where 
in Eq.(l.l) additional contributions have to be taken into account which are not 
products of the chosen set. But this extension is beyond the scope of our paper. Our 
principal interest is to show that besides of linear response, Mori's theory provides 
a powerful tool to find selfconsistent approximations for the time dependence of 
the Heisenberg operators G(t), thus leading to the possibility of calculating higher 
order response functions or expectation values for relaxation processes beyond linear 
response. 

Our paper is organized as follows: In section 2, after having given a short summary 
of Mori's theory, we present a formal expansion of the residual forces in terms of the 
dynamic correlation functions, and then discuss the approximation scheme for the 
points (i)-(iii). 

In section 3 we want to illustrate the formal results of section 2 from different points 
of view. First we will study the accuracy of the approximations by treating a simple 
model which can be solved exactly. Secondly, we will show how the formalism can be 
applied to a physical system, and will compare the results to other approaches. 



2 Expansion into powers of correlation functions 

2. 1 Summary of Mori 's theory 

To summarize Mori's theory we take the notation of [2]. The theory starts with a 
scalar product in Liouville space given by 



exp(iLt) (GG) = (exp(iLt) G) (exp(iLt) G) 



(1.2) 




(2.1) 



o 
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where B(t) denotes the time evolution in the Heisenberg picture. Projecting the 
observables G>(i) onto the space spanned by all the G v with 

PG v = {l-Q)G v = G v , (2.2) 



and onto the orthogonal space, one obtains the Langevin operator for the column 
vectors G(t) in the form 

t 

G(t)=iG(t)n - J G(t-t') {G\G)-\f\f{t'))dt' + f{t) , (2.3) 



where fit) is the column vector of the residual forces with the property 

Qf(t) = f(t) . (2.4) 

The abbreviations mean 

(G(t)% = E G -(*)^ ( 2 - 5 ) 

V 

((g|go- 1 (/i/(o)^ = E(g|g);a 1 (AI/m(*)) , (2-6) 

and the frequency matrix is given by the Liouville operator L 

n u , = J2(G\G)-t(G x \LG,) . (2.7) 

A 



The main point is that the dynamics of the residual forces f(t) are not governed by 
the Liouvillian, but by QLQ. 

The normalized correlation matrix S(i) 

3„„(f) = Y,(G\GU(G x \G,(t)) (2.8) 

A 



obeys the matrix equation 

t 

E(t) = -jz(t- t')(G\G)- l (f\f(t'))dt' . (2.9) 
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With help of the correlation matrix S(t), the Langevin equations (2.3) can formally 
be integrated to yield 

t 

G(t) = GE(t) + J f(t')E(t - t')dt' . (2.10) 
o 

It is this form which we will use in our following treatment. 

2.2 Expansion of f(t) in terms o/S(t) and static correlations 

To find an expansion for the residual force we first establish an exact nonlinear system 
of equations for /(f) in terms of the correlation functions (2.8). This will be achieved 
in two steps: We express f(t) by the time derivatives G(t), insert Eq.(l.l) which 
leads to products of G(t). Then the decomposition (2.10) of all G u (i) is used to get 
a closed system for f(t). 

Take the time derivative of Eq.(2.10), then using Eqs. (2.2) and (2.4) we find 

f(t) = QG(t)-f(t)®E(t) , (2.11) 

where we abbreviate the convolution of two functions a(t) and b(t) by 

t 

J a(t - t')b(t')dt' = a(t) b(t) . (2.12) 
o 

Substituting the derivatives (1.1) at time t 

[H, G,(t)} = £ V$G v {t) + E vl%G v {t)G x {t) + ■■■ (2.13) 

abbreviated by 

[H, G(t)} = V^G(t) + {G(t), G(t)} + ■■■ (2.14) 

into Eq.(2.11) we arrive at 

/(f) = iQV^ {G(t), G(t)} - if(t) ® (G\G)-\G\VW {G (t), G(t)}) + • • • ,(2.15) 
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where use has been made of the definition of E(t) (2.8), and of the relation 

Q(G(t) - /(f) <g> E(t)) = V^Q(GE(t)) = (2.16) 

which follows from Eqs.(2.10) and (2.2). 

Finally we insert the decomposition (2.10) of G(t) into Eq.(2.15). Then a closed set of 
equations arises which implicitely defines f(t) in terms E(t) and static correlations. 
To make it clear we write down the result for the case, where ■ ■ ■ vanish. 

One yields 

f (t)=iQvW {GE(t) + f(t) ® E(t), GE(t) + /(f) ® E(t)} (2.17) 
- if(t) ® (G\G)-\G\VW {GE(t) + f(t) ® E(t), GE(t) + f(t) ® E(t)}) . 

For simplicity we restrict to this case in the following. 

The implicit system (2.17) for f(t) is our basic result and will be the starting point 
for our approximations. To find an explicit expression for f(t) we will iteratively solve 
Eq.(2.17). Let us illustrate this procedure for a simple case. 

The expansion to be chosen depends on the magnitude of the matrix elements of 
V^ 2 \ Suppose that we can take all matrix elements of to be of the same order 
of magnitude of a smallness parameter e, 

V {2) ~ e . (2.18) 

Then we can write 

f(t) = ef^(t) + e 2 fW(t) + ... , (2.19) 

and comparing both sides of Eq.(2.17) we will find the explicit expressions for f^ n \t) 
in terms of correlation functions E(t) and products of operators G. 

The lowest order contribution yields 

f(t)=iQV^{GE(t),GE(t)} + --- , (2.20) 

where f(t) lies in the space spanned by G and the products GG, and the time- 
dependent coefficients are given by products of E(t). 

From a physical point of view the result (2.20) can be interpreted as a sum of coupled 
modes, whose time dependency is given by E(t). In Eq.(2.20) the projector Q acting 
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on the products GG projects out all contributions of the linear space spanned by 
the modes G v . Therefore the effect of the coupling of the S(i) can be very weak, 
although formally the bare interaction strength appears. This point will become 
more clear in our example in section 3, where we treat interacting spin waves. A 
further comment should be given as to the time scale for f(t) appearing in Eq.(2.20). 
In a macroscopic interacting system the matrix implies a summation over a very 
large number of products Z UfJj (t) E.\ a (t). So the phase factors of these products can 
produce a correlation time for f(t) which is entirely different from the relaxation 
times of Z(t). 

Iterating Eq.(2.17) one step further one obtains f^(t), which then results to be 
in the space spanned by G, GG, GGG. Although mathematically possible, we have 
found from the example discussed in section 3.1 that instead of calculating f^(t) it 
is more adequate to take a larger set of observables G which comprises the products 
GG, and use the lowest order approximation (2.20) to this set. Then one obtains 

f(t)=iQVW {GE(t),GE(t)} , (2.21) 

where f(t) is the column- vector of the residual forces arising for the set G. It includes 
the contribution of of the original set {G}. The relation between f{t) and f^(t) 
is discussed in appendix A. 

The simple expansion scheme starting with Eq.(2.18), or Eq.(2.20) respectively, will 
be illustrated in section 3.1, whereas an expansion of f(t) with 

V® = \/ (2) +eVf } (2.22) 
will be the basis for the treatment of our example in section 3.2. 

2.3 Expansion of G(t) in terms ofE(t) and static correlations 

Regarding the general decomposition (2.10) of G(t) into PG(t) and QG(t) it is clear, 
that inserting approximations for f(t) one obtains an expansion of G(t) into powers of 
S(t) and static correlations. For the case (2.18) the lowest order contribution follows 
to be 

G(t) = GZ(t)+iQV {2) {GZ(t),GZ(t)}(g)Z(t) + --- . (2.23) 

The main point of this result is, that for an interacting system, S(t) exhibits damping, 
so that the Fourier transforms of G(t) do not have divergent denominators. Further- 
more the projector Q again weakens the effect of the given interaction {G, G}, 
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as the part P{GG) is already contained in the dynamics of S(t). The structure of 
the result (2.23) is the same as it is known from an ordinary expansion of G(t) with 
respect to V<® . For [H {0 \G] = G, we would find Eq.(2.23) without Q, i.e. the 
full interaction, and S(t) being replaced by the transposal of exp [i V^t]. 

Once an approximation of S(i) is known, one can use Eq.(2.23) to calculate expecta- 
tion values for relaxation processes or higher order response functions. If one wants 
to improve the result (2.23), one can again take the larger set of observables G and 
use Eq.(2.23) for the corresponding quantities. 

2.4 Nonlinear equations for E(t) in terms of static correlations 

The correlation matrix S(t) is determined by Eq.(2.9), i.e. by the correlations (f\f(t)). 
Therefore use can be made of the expansion procedure for the residual force f(t) of 
section 2.2 leading to closed nonlinear equations for Taking the case (2.18) 

and substituting the lowest order contribution of f(t) (2.20) into Eq.(2.9) yields the 
following set of equations 

E(t) = iE(t)n 

t 

^ j^{t-t'){G\G)-\v^{G,G}\QV^{G'E{t'),G'E{t')})dt' . (2.24) 

o 

It is clear that higher order equations can be obtained, if higher order terms of f(t) 
are inserted into (f\f(t)). So we have found a systematic way to generate equations 
for correlation functions S(t). 

But we need not write down these higher order equations for S(t), because - as 
already has been discussed in section 2.2 - the higher order approximations to S(t) 
are included in Eq.(2.24), if used for an extended set of observables G = {G, GG — 
(GG) , GGG — (GGG), •••.}. with new parameters and and the corresponding 
correlation matrix 

Inspecting the result, (2.24) shows that it has the structure known from mode-mode 
coupling approximations. In our case we have obtained equations - as mentioned in 
the introduction - with bare interactions [5-7]. 

As a check of the validity of the approximation (2.24) one can take the time depen- 
dence G(t) resulting from (2.23) with the approximate S and calculate the correlation 

A(t) = (GIG)- 1 (G(t)\G(t)) (2.25) 
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which for the exact dynamics G(t) is one. In the discussion of the example given in 
section 3.1 we have found that the deviation of A from 1 is directly related to the 
accuracy of the approximate S. 



2.5 Static correlations 



The equations of motion for S(t) Eq.(2.24) contain static correlations. They can 
be viewed as given parameters and can be taken from any static theory which is 
available. In this section we want to point out that in many cases however, one can 
close the dynamic equations (2.24) by a set of equations for the static correlations 
which appear in S(t), so that all quantities - at least in principle - can be determined 
in a selfconsistent way. The key point is that the approximations obtained for the 
observables G(t) in the Heisenberg picture are expressed in terms of c-number time 
functions and products of operators G at time t — 0. So the well known KMS- 
condition for the Heisenberg operators G at imaginary times if3 

(GA) p = (AG(0))p , (2.26) 



which can be cast into the relations 



(A\G)= J — {[G(u),A^]) p (2.27) 

— oo 
oo 

{AG) P = J duj(e^-l)- l {[G{u),A}) p , (2.28) 



directly connects static and dynamic correlations. Here G(u) denotes 

oo 

G(u) = — [ G(t)e ZUJt dt , (2.29) 

27T J 



and it has been assumed that the observables G are chosen to be orthogonal to 1 and 
no other constants of motion are to be projected out. Therefore we must only show 
that given an approximation for G(t), or f(t) respectively, Eq.(2.27) and Eq.(2.28) 
will provide us with a closed set of equations for the relevant static correlations 
appearing in H(£). To this end we will suppose that the operators G obey a Lie- 
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algebra P"P1 



[G, G] = X] • • • G + CI 



(2.30) 



leading to 



[G,GG] = J2 ■■■ GG + J2-G (2.31) 
[GG,GG] = J2 ■■■ GGG + J2 - GG . (2.32) 

For clarity in the presentation we refrain from writing down the structure constants 
explicitly, but focus on the observables. Let us first discuss the expansion scheme 
y( 2 ) ~ e of section 2.2 with the lowest order approximation for G(t) found to be 

G(t) = GE + iQV {2) {GE,GE}®E . (2.33) 

The corresponding equation for H is given by (2.24). If we apply the Kubo identity 

P(G\LA) p = ([Gft ) A])p (2.34) 



to its frequency and memory matrices and carry out the commutators with help of 
Eqs.(2.30) and (2.31), and (G)p = 0, then we see that E(t) depends on the static 
correlations (G\G), (GG), (G\GG) only 

E = E{t,{G\G),(GG),{G\GG)) . (2.35) 



To connect these static quantities we start with two relations which are independent 
of the approximations. Regarding LG = V^'G+V^ {G, G} and calculating (G\LG) 
from the Kubo identity with Eq.(2.30) we arrive at 

C = f3{G\V {1) G)+f3{G\V (2) {G,G}) . (2.36) 



On the other hand we calculate (G^G) applying Eq.(2.28) to A = G^ . Inserting (2.33) 
and the equation of motion for E(u) we arrive at 

(GtG) = (G\G) J du-gg^Eiu) . (2.37) 



2 The matrix C can originate from the fact that the G have been chosen with the property 
(l\G)f, = {0)0 = 0. 

3 If non-Hermitian operators G are used, the operators G^ are assumed to be elements of 
the space spanned by the set G. 
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Two further relations are obtained, if we explicitly use the approximation (2.33) 
for G(t) or G(u) respectively. We regard Eq.(2.27) for A = {GG)^ and insert the 
Fouriertransform of (2.33)0. Then respecting the commutators (2.31) and (2.32) we 
find equations of the following structure 

(G\GG) = E ' ■ • ( GG ) + E ' • • (GGG) + ■ ■ (GG){G\G)-\G\GG) , (2.38) 

where we have omitted the coefficients which are given by the structure constants of 
the Lie algebra and matrix elements of 




(2.39) 



To close this system of equations with (GGG) we regard Eq.(2.28) for A = GG 
yielding 

(GGG) = ■ ■ ■ (GG) + ■ ■ ■ (GGG) + £ • ■ • (GG)(G\G)-\G\GG) (2.40) 

with coefficients analogous to the coefficients appearing in Eq.(2.38). The only dif- 
ference is, that the denominators j3u in Eq.(2.39) are replaced by {er w — 1). 

These four equations (2.36), (2.37), (2.38), (2.40) for (G\G), (GG), (G\GG), (GGG) 
now in principle allow for a determination of the static correlations in terms of func- 
tionals of the dynamic correlations S(a>), where it is clear, that the solvability must 
be discussed in each application. Although these equations are extremely compli- 
cated, one can think of an iterative solution starting with suitable static values in 
the dynamics of H. Another simplification occurs, if the temperature is high enough 
so that tanh flu can be replaced by j3u. Then it holds 

(G\G) = \([G\G] + ) (2.41) 
(G\GG) = 1 -([G\GG] + ) , (2.42) 

and the total system reduces to a set for the correlations (GG) and (GGG). 

We emphasize that the basis for our discussions on the statics was the result (2.33) 
for the Heisenberg dynamics. If we now consider higher order contributions to G(t), 
we have found two possibilities, as was pointed out in section 2.2 and 2.3. But as 
to the static correlations these two treatments are not equivalent. If we improve 

4 The correlations (G\G) cannot be calculated from Eq.(2.27) with A = G, as in this case 
Eq.(2.27) identically holds, if the equation of motion for E(a>) is inserted. 
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Eq.(2.33) by a residual force f {2) {t), so that G(t) is an element of the space spanned 
by G, GG, GGG, then it is not difficult to see, that the static equations resulting 
from (2.27) and (2.28) are no longer closed. However, if we use the alternative and 
extend the set of observables to yield a set G which includes the products GG, then 
the set of equations for the static correlations appearing in S(t) again can be closed 
(appendix A). This also indicates that the second way is more adequate. 

So far we have discussed the equations for the static correlations which result from 
the expansion scheme for the residual forces with ~ e. In general the static 
relations to be obtained will depend on the iteration procedure for f(t) which will 
be chosen. In our example of section 3.2, where we treat the lowest order of a case 
V® = V (2) +eVf ) , we will see that a closed set of equations can be derived again. 

The considered approaches leading to closed sets for the static correlations have the 
common feature that they are valid to the same degree of accuracy as the dynamics 
in the Heisenberg picture is correct. In this sense dynamics and statics are treated 
at the same level of approximation. 



3 Illustrations and discussion 

In this section we will illustrate the general formalism of section 2 applying it to 
two examples. First we will study general features of the approximation scheme. To 
this end we choose a model as simple as possible: We take a spin system with long 
range interactions in the high temperature limit, so that we need not handle with 
the problem of the static correlations, and can obtain both, the exact solution to the 
dynamics and the approximations in analytic form. As a second example we will treat 
a Heisenberg ferromagnet at low temperatures. We will show that one can find an 
expansion of the residual forces in terms of spin operators which leads to meaningful 
results for the Heisenberg dynamics and the correlation functions, so that we can 
make contact to other approaches. 

3.1 Exactly solvable model in the high temperature limit 

Our model to be considered is a Heisenberg spin system with long range interactions 
in the high temperature limit, where we can carry through the general expansions 
of sections 2.2 and 2.3 to high order explicitly. We will find that the lowest order 
approximation to the residual forces leads to a coupling of two modes only, where 
one of them is a trivial constant of the motion. So the system might be a very special 
one, but we think that it correctly gives insight into the problems of time scales and 
convergence. 
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The Hamiltonian of our system reads 



w = _ 7^p- % = -^- 4, (3 ' 1) 



where 



N 



S q = Y,e im % (3.2) 



i=i 



S — S q= o , (3.3) 

and the N spins (s = are located at lattice sites In the high temperature limit 
the scalar product reduces to 

(A\B) = (Trl^Tr^S), (3.4) 



so that static spin correlations can be evaluated. 

— * 

In choosing a set of observables {G} we take S q , but substract the components parallel 
to the constant of the motion So, so that the dynamic correlation functions E(t) will 
decay to zero (c.f. appendix B). Therefore we take the set of observables to be 



— * — * 

Go = Sq 

G q = S q -^(S -S r 1 ((S ■ S q )S + S (S ■ S q )) q^O (3.5) 
with So • S q — Y, a SoS q . Then it holds 

LG q = i -L(G q xG -G xG q ) q^O , (3.6) 

which means, that the conditions of section 2.2 

LG = V {1) G + V {2) {G,G} (3.7) 

are fulfilled, where vanishes and just couples G q to G . From the symmetry 
of the Hamiltonian it directly follows that the correlation matrix S(t) is diagonal with 
respect to q, q' and cartesian components a, a' = x, y, z, so that the decompositions 
(2.10) specialize to 
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q = : G (t) = G ■ 1 + (3.8) 
q^O:G q (t) = G q E(t) + f q (t)®E(t) , (3.9) 

where 

E(t) = (GJIGJ)- 1 (G£|G£(t)) (3.10) 

for g 7^ does not depend on g and a. So we are ready to study the expansion of the 
residual forces f q (t) with an expansion parameter e ~ J. 



3.1.1 The first order approximation 

We now take our model to carry through the expansions of section 2 for the case 
yCO = 0, V {i) ~ e ~ J iterating the Eq.(2.17) for f q (t). Although J will be a formal 
expansion parameter the accuracy of the approximations will not be determined by 
the smallness of J but rather by a condition on time Jt <C 5, as for = we just 
have one time scale given by J -1 . Nevertheless in a high order approximation the 
bound S may be so large that the "short time expansion" covers the whole range of 
physical interest. We will study this problem of time range in detail, calculating the 
correlation function S(t) to general order. 

— # 

The lowest order approximation for the residual force f q (t) follows from (2.20), (3.8) 
and (3.9) to yieldF] 

f q (t) = -^=(G xG q -G q x G )~ = iLG q E , (3.11) 

where just two modes E q = H and S g=0 = 1 couple, so that the force f q (t) for all 
times is proportional to one fixed element in the space of the GG. 

The dynamics of the Heisenberg operators G q (t) corresponding to the approximation 
(3.11) are obtained from (2.23) to give 



G q (t) = G g E(t) + -L(G x G q - G q x G )E(t) ® E(t) 



= G q E + iLG q E ® 5 . (3.12) 

— * 

In this expansion the coefficient E(t) of G q is exact, whereas the coefficient E ®E 
of iLG q will be modified by higher order terms, if these are projected onto LG q . 



As (G q \LG q ) = the projection operator Q can be ommitted. 
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The equation of motion for the correlation function E(t) is given by Eq.(2.24) and 
simplifies to 

Z = -m 2 Z®Z , (3.13) 



where 



is independent of q and a. The solution of (3.13) for S(i) can be expressed by the 
Bessel function J\. Taking the Laplace transform of (3.13) one finds 

*')(.) = 2 Z£±^±M , (3.15) 



or 



S(i)(t) = 2^^, r = 2v^t (3.16) 



respectively, where the index denotes the order of the approximation. This result for 
S^(i) is compared to the exact solution in Fig.l. One sees that S^(i) is close to 
the exact correlation function, as long as Jt is smaller than the zero of S(t), i.e. as 
long as it holds 

Jt < 1 . (3.17) 



At the end of section 3.1 we will show that the range of validity (3.17) for E^(t) can 
be estimated from E^\t) and (3.12) alone, without knowledge of the exact solution. 



3.1.2 Higher order approximations 

As there are more possibilities to derive higher approximations we will first study 
the second order approximation in some detail and then give the general results. 

Iterating Eq. (2.17) for f q {t) up to order J 2 we find for the residual force 
J 



f q (t) = -^=(G xG q -G q xG )~ 
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Fig. 1. Approximations H^'^r) compared to the exact correlation function H(r); a: n = 1, 2, 
exact; b: n = 3, 4, exact; c: n = 5, 6, exact; d: n = 7, 8, exact. 

+ ^ Q {G x (Go xG r G,x G ) - (g xG q -G q x G ) xG„}»®» 
= iLG q E + g(iL) 2 G 9 S ® S , (3.18) 

and for the Heisenberg dynamics of G q {t) 

G q (t) = G q E + iLG q E(g)E + Q(iL) 2 G q E(g)E(g)E . (3.19) 



The dynamics of E(t), however, keeps to be E^(t), as the new added term in f q (t) 
is orthogonal to f q (0) = iLG q , 

(f«\Q(iL) 2 G a q ) = , (3.20) 



so that the memory function does not change. To improve E^\t), at least third 
order terms in f q (t) have to be considered which have a non-vanishing projection 
onto iLG^. 

An alternative approach to the second order result described above can be given, if 
we keep the basis we have obtained for f q (t) and G q (t) in (3.19), but use a formulation 

— * 

such that the factor of iLG q will not be changed by higher order terms. This means 
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that higher order terms must be constructed to be orthogonal to the products iLG q = 
J/VN(G x G q — G q x G ) (c.f. appendix A). 

For this alternative approach to the second order approximation we extend the set 
of observables G to the larger set {G} , which comprises the occuring basis vectors 
of products GG, i.e. iLG q : 

{G} = {Gq, G g> i = G q , G Q) 2 = iLG q }. (3.21) 



The corresponding projection operator onto this space is denoted by P = 1 — Q. 
Then the general formalism of section 2 applies to the set {G}, as it holds 



iLGq t x — Gg j2 

iLGq, 2 = ^=(G xGq, 2 -Gq, 2 xG ) . (3.22) 



N 

Thus the condition (1.1) with V {1) ^ and V (2) ~ J again are fulfilled, and the 
lowest order approximation (2.20) for the residual forces j/j = {/i^i,/^} reads 

£,i(*)=0 (3.23) 
/; 2 (i) = 4=Q(GoxG,, 2 -G 5 , 2 xG )H 22 , (3.24) 



where the matrix of the correlation functions is defined by 



v,fi = l,2 . (3.25) 



(G a \G a ) 

Hence the Heisenberg dynamics for G q {t) = G q> i(t) follow to be 

Gq(t) = G qt i{t) — Gq^Hn + G gi 2^21 + fq,2 ® "21 

= GqE u +tLG q E 21 + Q(tL) 2 Gq(E 22 ®E 21 ) . (3.26) 

One sees that the result (3.26) for G q (t) is spanned by the same basis vectors as it 
was in Eq.(3.19). Just the time-dependent coefficients have changed.^ 



6 If one extracts the residual force f g (t) with respect to {G} from the result (3.26), one 
finds that f q (t) also has the same basis vectors as the approximation (3.18). It is not difficult 
to see, that the memory functions calculated from this f q (t) just lead to the exact relation 
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In order to determine we now write down the equations of motion for the matrix 
( c -f- Eq.(3.22)) which follow from the frequency matrix (2.7) and the memory 
matrix given by {f q2 l/^O)- We obtain 



—11 = —12 

S12 = -c 2 Sn - C3H12 ® H 22 (3.27) 

— 21 = — 22 (3.28) 

S22 = -C2H21 - c 3 (S 22 <8> S 22 ), (3.29) 

where C2 = m 2 and C3 = (m 4 -mj)^, and m 4 denotes the fourth moment. 

This system can be solved by Laplace transform to yield for = S n 



h« 2 )(,) = ^ v — — *_l . (3 30) 

s 

c 2 

The corresponding result in time S^-'(t) is shown in Fig.l. As compared to S^(i) 
the range of validity of S^ 2 ) (t) has obviously increased. 

The alternative approach for the second order approximation can be extended to 

— * — * 

general order. Instead of iterating the Eq.(2.17) for f q (t) and calculating f q (t) up to 

— — * — * — * 

order J™, we choose a set {G} given by G q ,iLG q , ■ ■ ■ (iL) n G q , and take the lowest 
order approximation for the residual forces. The details are given in appendix B. Here 
we just list the results. To write them down, it is expedient to use an orthogonal basis 
G q)U in the space |g| introduced by 



Gg,l - G q 

G q ,2 = iLGg : \ 



I Via \C a s 



G q ,v = iLG q ,„- 1 + —% — * rG q , v -2 z/ = 3,---,n . (3.31) 

V^q,v-2 \^q,u-2 ) 

Then the residual forces read 



/,,„(«) = v = l,---n-l 

fq,n{t) = G^n+l— nn , (3.32) 



' (G°|G?) ® — 
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where G q , n +i is defined by (3.31) with v — n + 1. The Heisenberg dynamics result to 
be 



G q (t) — G qt i(t) — Gg^Z.vl + G q , n+ i^ 



nn <y "nl 



(3.33) 



i/=i 



The equations for the correlation matrix 



(3.34) 



are found from the frequency and memory matrix. The explicit solution for the 
Laplace transforms of E ul and S nn are 



=■1/1 (s) = (-1) i/ = 2,---,n , 

C l " " ' Cj, 



where = Hn is given by 



B r , 



Bn-1 + \ V-Bn-1 



S (n) (s) 



-Br; 



+ 4c 



n+l 



2 (-iy 



1 c n+l D 2 



Cl • • • C r , 



(3.35) 
(3.36) 



(3.37) 



and the A v (s), B v (s) are polynominals in s defined by 

A v = s + C v A u - 2 A = A x = 1 
B u = s B v _i + CV -Bo — 1 Bi — s 



v = 2,---,n (3.38) 



with 



ci = l 



(^g,i/-ll^g,i/-l) 



u = 2,---n, a = x, or y,z . (3.39) 



The functions 'E^(s) for complex s are holomorphic for Res > and have the correct 
property 

ReH (n) (s) > 0, Res^O . (3.40) 
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0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 



CO CO 

Fig. 2. Spectral densities of the approximations ReH( n )(iu), compared to the exact density 
for n = 2, 3, 12, 13. An averaged spectrum is shown for n = 12. Frequency and spectrum in 
units of (4m2) 1/ ' 2 and (4m2) -1 / 2 respectively. 



The results for 5^(t) and the exact solution H(t) are plotted in Fig.l. One sees that 
the time region, where ^ n \t) is a good approximation to increases with n. The 
spectral densities 1/2 dtE(t)e~ lult = ReS^(s = ia;) are shown in Fig. 20. The 
intervals where ReS( n )(iu;) identically vanishes are due to the square root occuring 
in S( n )(s) and are a consequence of the coupling of just two modes. Therefore in the 
case of many modes one can expect that the values ReS^ n ^(io;) ^ are no longer 
restricted to finite intervals which cause the asymptotic oscillations for long times in 

At the end of this section we will discuss the validity of the approximations from 
a different point of view. Our approach does not only lead to equations for the 
correlation functions. The basis was the expansion of the residual force and the 
resulting Heisenberg dynamics G q (t). From this we will be able to write down a 
necessary condition for the accuracy of the n-th order approximation for and 



7 If one wants to get approximations which are valid for all times ^ Jt ^ oo one 
can average the spectral densities ReH( ra )(zu;) with a width 5 which smoothens the intervals 
ReH(")(icj) = 0. An example for n = 12, 1/tt f duj'ReE( n \iuj')5/[{uj-uj') 2 +5 2 } = ReH( n )(cH 
iu) is given in Fig. 2. 
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G q (t). Let us introduce the correlation at equal times 



A t = v y'\*: 3.4i 



— # — # 

which for the exact G ? (i) must be one. Inserting G q (t) from (3.33) we find 

n 

(t) = E Si„(*)S„i(t) + Cl c 2 • • • • c n+1 (S n „ <g> ~ nl ) 2 (3.42) 

!/ = l 



where the orthogonality of the and the definitions (3.34) and (3.39) have been 
used. The expresion (3.42) for A*™) can be simplified further with help of the equations 
of motion for E UfM (t), or the explicit solutions (3.35) and (3.36), e.g. 



+ — 

c 2 



(3.43) 



If the n-th order approximation is appropriate we must have 

A (n) (t) ~ 1 (3.44) 



which gives a restriction to the region of time, as for t = it holds A( n \t = 0) = 1. 
The deviation of A( n )(t) from 1 in Eq.(3.42) comes from the finite basis in Liouville 
space entering into G q (t) and the approximation for the correlation functions, which 
is related to this subspace. Thus A^ n \t) ~ 1 requires that the restricted Liouville 
space in G q (t) as well as the approximate E^ n \t) are sufficient. In Fig.3 we have 
plotted A^(t). Comparing the deviations S^(t) from E(t) to A^ n \t) one observes 
that E^ n \t) is a good approximation to E(t) 

sW(0~S(t) , (3.45) 



as long as it holds 



A (n) (t)>^ . (3.46) 



From this we conclude that one can take the condition (3.44) as a measure for the 
accuracy of the approximations without knowing the exact solution H(t). 
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3 6 9 12 15 

X 

Fig. 3. Accuracy of the approximations; norm A^ n \r) compared to the deviations of 
H("0(t) from the exact correlation function S(r) for several orders of the approximation 
n = 1,2,... ,13. 

3.2 Heisenberg ferromagnet at low temperatures 

Introducing the Heisenberg ferromagnet we want to show how our general method 
can be applied to a system where temperature dependent static correlations occur. 
We will not present more refined results, but rather want to study the selfconsistent 
method by comparing the results to those of other approaches. To this end we will 
simplify our equations for the dynamic and static correlation functions, so that we 
can recover standard results for temperature dependent spin wave frequencies and 
dampings. In applying our formalism we keep the spin operators and do not transform 
to Bose operators. Thus the Heisenberg dynamics of the spin operators are expressed 
in a basis of spin operators at t — which may have some advantage for further use. 

The Hamiltonian of our system is given by 

i i^j 

= -HS Z - l/N Y,J q S q - H -> +0, J i3 > (3.47) 

— * — * 

where N spins S are located at sites Ri in a lattice and S q is given by Eq.(3.2). We 
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want to treat this system for low temperatures. To this end we introduce a parameter 
e to be the inverse of the spontaneous magnetization 



<f>) =v-\T) , (3.48) 



and scale the time to 

r = at , (3.49) 

and the spin operators according to 



Sf = -£= • (3-50) 



The operators 5S* — S* — (S*) will not be affected. Then the equations of motion 
for the scaled spins read 

±8Sl = ~ - J2)S 2 SU (3.51) 

±Sf = T2t(J o -Ji)Sf±e^Y,( J i-2-J2)SS0t2 , (3-52) 

where we have used the shorthand notations q\ — > 1, q\ — q 2 — > 1 — 2. The parameter 
e(T) will not approach zero, except for the case S — > oo, but one can directly see from 
(3.51) and (3.52) that the zeroth order contribution in 8S((r) and S^(r) gives the 
correct spin wave approximation in lowest order [8]Q, so that introducing e will allow 
for a low temperature approximation. This will be confirmed by the final results. 

We want to apply the theory of section 2 to the dynamics with respect to r, which 
means that we formally have a Liouvillian L = L/a, or Hamiltonian TC = T~C/<J, 
respectively. As a consequence we have a scaled temperature $ with = (3H.. 
Choosing the set of observables G to be all modes 5S* and Sf- 

{G} = {6S*Sf} , (3.53) 



and inspecting the equations of motion (3.51) and (3.52) for G one sees that it holds 

LG=V {1) G + V {2) GG (3.54) 



8 As referred to time t we have the renormalized spin wave frequencies oj\(T) = 2<r(T)(Jo - 
Ji). 
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with 



= \/ (2) + eV} 2) 



(3.55) 



so that our condition for an expansion of the residual forces is fulfilled. Furthermore 
the G, 1 form the Lie-algebra 



St,S; =N5 1: - 2 + e5S z 1+2 



=C °l+2 



(3.56) 
(3.57) 



which will be important for the statics. 



3.2.1 Residual forces and Heisenberg dynamics 

We start with the exact equations (2.17) for the residual forces which follow for our 
special dynamics. Regarding that the static and dynamic correlation matrices are 
diagonal with respect to wave numbers and components a = z, +, — , and introducing 



5i(r) 



Sf(r) 



(SSf 


5S*(t)) 


(5Sf 


8Sf) 



(S?\St(r)) (Sf|Sf(r)) 



(s?\sf) 



(st\st) 



(3.58) 
(3.59) 



we have 9 



ffc) = -| E( J i-2 - Ji)Q (4~s 2 - + / 2 - ® sj) (s+ 2 s+ 2 + /+ 2 ® s+_ 2 ) 



iV 



+/^lE (Jl " a " Ja 



iV^ («5Sf|«5Sf) 

(4-s 2 - + / 2 - ® s-) (5+ 2 s+_ 2 + /+ 2 ® ~+_ 2 )) , 



(3.60) 



and 



/±(r) = ±e| £( J!_ 2 - J 2 )Q (*S«S* + /| ® S£ (Sf_ 2 Hf_ 2 + /± 2 ® S±_ 2 ) 



? ± , 2i ^ (Ji_ 2 - J 2 ) 



9 The convolutions refer to r. 
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•(^|(^S5 + ^®^)(^_ 2 Sf_ 2 + /f_ 2 ®Sf_ 2 )) . (3.61) 

One directly sees that the transversal forces fi (r) start with a first order contribution 
in e whereas the longitudinal forces /-f (r) have a term of zeroth order. This structure 
of Eqs.(3.60) and (3.61) allows for a solution in powers of e. Restricting to the non- 
vanishing lowest order approximation 10 we find 



= -| E S£ 2 (Ji_ 2 - J 2 )0i 2 (S + , S-) (3.62) 



/^(r) = ±e| £(^-2 - J 2 )Q {*5^_ 2 ^Sf_ 2 



2i 



N 3 



£ Q (4"4 + - 3 ) S± 2 ( J 2 _ 3 - J 3 ) (023 ® 5 2 ) ~f_ 2 , (3.63) 



where </>i 2 is a functional of the correlations H + , S and can be expressed by Laplace 
transform to yield 



<M S ) = tt^At r • ( 3 - 64 ) 









1(a) 


2t (5Si\SaSt 3 ) 
(SSf\SSf) 


-(J^s- J3)(S 3 -S+_ 3 )( S ) 



The Heisenberg dynamics of the spins follow from the general relations (2.10). Re- 
garding the symmetry of the correlation matrices we find 



SSfa) = 5SI~{{t) + /*(r) ® S*(r) (3.65) 
5 1 ± (r) = 5 1 ± S^(r) + A ± (r)®S±(r) . (3.66) 

Inserting the approximations for the residual forces (3.62) and (3.63) into (3.65) 
and (3.66) respectively, yields the corresponding approximations for the Heisenberg 
dynamics. One sees that the longitudinal components SS^(r) move in a subspace 
spanned by SS^ and all products S 2 Si_ 2 , whereas the transverse components are 
restricted to a space spanned by Sf and all the products <5S'|S'j b _ 2 and S^S^sS^ 
respectively 



5St(r) = • • • <^ + » • S 2 St 2 (3.67) 

2 



io p or j± ( r ) ^he g rs ^ orc [ er contribution is necessary as to obtain spin wave damping, 
whereas the first order correction to /f (r) would just modify Sf (r) quantitatively. 
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Sf(r) = ■ ■ ■ Sf + e£ • • • 6S* 2 S?_ 2 + e £ • • • 

2 2,3 



(3.68) 



The time dependent coefficients are given by functionals of the dynamic correlations 
S z (r), S ± (r), and the exchange parameters, and static correlations. It should be 
noted that the results for SS z (r) and 5 ,± (r) have the same subspaces as would have 
been obtained by a simple perturbational expansion of the dynamic equations (3.51) 
and (3.52). The essential difference of (3.67) and (3.68) to such an expansion is, that 
the coefficients of 5Sf and Sf in (3.67) and (3.68) comprise all contributions of the 
simple perturbational series which have a projection onto SS*^^. As compared to 
a Holstein Primakoff approach^], 5 ,± (r) corresponds to a perturbational treatment 
of the four magnon interaction, if in addition, the Bose operators are summed up to 
give the spin operators, and the coefficients will be renormalized. 



3.2.2 Dynamic and static correlations 

The dynamics of the correlation matrix 5 is governed by the Mori equations (2.9). 
Regarding the symmetry of the Hamiltonian we get 



i. S | = _ 7 ll 0S * (3.69) 
Asf = =FiwiSf-7r®E r , (3-70) 

GST 



where we have introduced the memory functions 



^ w =$S^ = WM) <fcJlf(T)1> ' (3 ' 72) 

and the frequency 

Wi(T) = N (PiS+lSt))' 1 . (3.73) 



11 Expanding the coefficients in (3.68) into powers of e with e(T = 0) = 1/S and inserting 
the expansions of into Bose operators a, a) into (3.68), then up to order 5 _1 the result 
for Sf(r) coincides with the dynamics of a(r) calculated with the four magnon interaction 
as a perturbation. 
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In the following we will use the commutator forms for 7II and 7 J ~[ T7 1. Inserting the 
approximations for /f(r) and /1 (r) (3.62), (3.63) into (3.71) and (3.72), we obtain 
the following expressions for the memory functions: 



7) (r) = |f E( J i-2 " J2) A 12 (T)<j> 12 [~+, S~] (3.74) 
iV 2 

iV 2 

+ e (f ) 2 E( J i-2 - ^)(J 2 - 3 - h)C 12 ,z{T) (0 23 ® H2) S+_ 2 , (3.75) 

2,3 

where the temperature dependent coefficients are given by 



^( T ) = ^4tt^ - 0M± a >) (3-76) 



B * m = msn i {M - e ™ - i iBr ) (377) 

r< (r^ — (f^fljg ^2-3) R /^x 



+ mW){ iA2 - i+s "' )N{S;s - 3) wist) 

+ £ <S 3 -iS|_ 3 S+ 2 > + £<5 3 -S+ 3*SJ>1 . (3.78) 

Inspecting the time dependence of the memory function, one sees that 7" and 'j ± are 
determined by S + , H - , or S + , H - , H z , respectively. Thus we have a nonlinear coupled 
set of equations for This set of equations is very complicated, as it is 

nonlinear and contains static correlations which must be known, or determined from 
the dynamic equations. For this reason we first state, how the static correlations in 
principle can be obtained in a selfconsistent way. Then we calculate their leading 
terms for T — > and show, how the equations for S z and S + reduce to the standard 
forms. 

The static correlations entering into S z and are found from Eqs. (3.76)-(3.78) to 
read 

{SS'\SS"), (6S Z 6S Z ), (S+\S + ), (S-S+), 
{5S Z \S-S + ), {S + \5S Z S + ), (S + \S-S+S + ), (S-S + 5S Z )' { ' ' 



12 We avoid Mori products of the form (GG\A). 
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If additionally one takes into account the correlations 

(5S Z 5S Z ), (S-S-S+S+), (S-S + 5S Z 5S Z ), (S~ S~ S + S + 5S Z ) , (3.80) 

then one can find a closed set of equations for this extended set (3.79), (3.80). The 
procedure is similar to the considerations in section 2.5. The difference is, that our 
iteration with + eV^ has lead to residual forces which also have contributions 
of the form GGG, so that the commutators are modified. The details are given in 
appendix C. 

For our purpose, however, it is sufficient to have the static correlations for T — > 0. 
To this end we iterate the system of static equations with respect to e, as the first 
order contributions result to be smaller than those of zeroth order. In zeroth order 
the coupled system of dynamic and static correlations for Sf , Sf and the correlations 
(3.79) can be solved exactly. Then the expressions of first order follow in a straight 
forward manner. As those terms we need are exact, they can also be obtained in a 
different way. Introducing 

a>i = 2(J - Ji) (3.81) 
m = -=— (3.82) 

we explicitly find for the transverse correlations 13 

(SlJt) = Nm + em(l + m)/3^(aJ 2 - Cj^ 2 )n 2 + C(e 2 ) (3.83) 

2 



{st\st 



-ni-2 



-e jrai_ 2 (l +ni_ 2 )-^2(a)3 -wi_ 3 )n 3 - ^^M 1 +^2-3) 

+ ~K~T S[ 7l s( 1 + ni ~ 2 + n2 - 3 ) ~ n l-2^2-3] 



N 3 

0J 2 — £b\-2 + 0)1+3-2 — £>2-3 
UJ 3 — U)\- 2 — & 2 - 3 , 

(S\ \S 3 S 2 _ 3 Si_ 2 ) , , 

— = (s ' 2 - 1 M 3 

+ e \ (£3,2-1 + S 2 , )n 3 (l + n 3 )f3 



+ 0(e 2 ) (3.84) 



13 The results are given in zeroth order, and up to first order, where they will be needed. 
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• ~~ ^3-4)^4 - n 3 (l + ni_ 2 + n 2 -3) 

4 

+ [n 3 (l + ni_ 2 + n 2 _ 3 ) - ni_ 2 ^2-3] 



^2 — ^1-2 + ^1+3-2 — ^2-3 



+ C(e 2 ), (3.85) 



CJ3 — — ^2-3 

whereas the longitudinal correlations are calculated to yield 

(5S z _ 1 5S z 1 ) = J2n 2 (l + n 1 ^ 2 ) + 0(e) (3.86) 
2 

{Sl 2 StiSSZ) = -Nn 2 (l + m_ 2 ) + 0(e) (3.87) 

1*5?) = £ £ " : + ^(e) (3-88) 

= ~N £ - ""'"I : + 0(6) • (3.89) 

Pl^l-2 - ^2) 

All the static correlations (3.83)-(3.89) listed above still depend on the spontaneous 
magnetization cr(T) = e -1 . This can be determined from the condition X^s* ■ s« = 
NS(S + 1) which can be cast into the form 

= (a-S) + (v-S)(j^ E<^A + >(^) + 25+l) 

+ ^2 E (<^ 1 ^>(a) + 25(51^) ((7)) (3.90) 



and leads to the magnon result 

= .<? - 



Once the statics are known we now can go back to the dynamics of S 2 and S ± , 
calculating the lowest order contribution to the memory functions (3.74) and (3.75). 
For 7^ (s) we find the zeroth order result 



ni-2 - n 2 



2 (U) 2 ~ U1-2) (S ~ i{£0 2 - 0)1-2)), 

(£2 - tDi_ 2 )(wi-2 - n 2 ) 

2 S - i(0> 2 - £1-2) 



E 



(3.92) 
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which leads to 



whereas the memory function 7j_(t) is of second order and reads 



7i (r) = (— ) M }Jw 2 - wi_ 2 ) = = = 

J\ ^ UJ 3 — LU 2 -3 — LO1-2 

•[n 1 _ 2 n 2 _ 3 (l + n 3 ) - (1 + m_ 2 )(l + n 2 _ 3 )n 3 ]e^ 3 "" 2 - 3 ^ 1 - 2 ) T . (3.94) 
The spin wave frequency (3.73) up to second order yields 



wi(T) =&! - ^- X)(^2 - ^1-2)^2 
^ 2 

~~ (t?) 2 XX^ 2 - £1-2) ^n 3 [n 2 _ 3 + n 2 (l + n 2 )/3(u; 3 - u>i_ 3 )] 
iV 2 3 

tf (T = > . (3.95) 



/K>i 

Comparing our results for T — > to those of spin wave theory and regarding 



™ 1 = t-2a(T)(J - Ji) (3.96) 
Pu 1 = P-2<r(T)(J -J 1 ) (3.97) 

we see that S 2 (t) coincides with the correlation function obtained by Lovesey [8], 
provided we take our temperature dependent spin wave frequency 2<t(T)(J — J±) at 
T = 0. The memory function Ji(t) = 7^ (r) (dr / dt) 2 (3.94) agrees with the Holstein 
Primakoff result for the four magnon interaction [5], if again cr(T) — > S is used. 



It is clear that the presented procedure for the Heisenberg ferromagnet is rather 
complicated. Just for calculating spin wave frequencies and dampings one would 
choose the usual way. But we think it is important to have shown that one is not 
forced to expand spin operators into Bose operators, as to obtain the correlation 
functions at low temperatures. In this context our main result is, that the Heisenberg 
operators Si(t) (c.f. Eqs.(3.65) and (3.66)) can be expressed in terms of correlation 
functions and spin operators at t — 0. This point will be useful, if one wants to derive 
equations of motion for expectation values of spins which go beyond linear response, 
as in such a case the expectation value of (Sg)(t) cannot simply be replaced by the 
expectation value of one Bose operator (a q )(t) 
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4 Conclusion 



It has been shown, that incorporating the time evolution of products of observables 
(GG)(t) = G(t)G(t) into the frame work of Mori's theory, one can give selfconsistent 
approximations for the residual forces. These can be expressed in terms of time- 
dependent correlation functions and operators at t — 0. In this way it is possible 
to deduce approximations for the Heisenberg dynamics of the observables as well 
as for the dynamic and static correlation functions. We have tested this approach 
comparing the approximations to the exact solution of a model and to the theory 
of interacting spin waves. From this comparison we conclude that the selfconsistent 
approximations in Mori's theory can be successfully used to attack the dynamics of 
correlation functions and Heisenberg operators of interacting systems. 

The dynamic equations for the correlation functions we have found, are of the type 
of mode-mode coupling equations, but the coefficients involved are given by bare 
interactions. In addition, our approach allows for a quantitative estimation of the 
validity of the presented mode-mode coupling approximation. A point of further 
investigation would be, how the expansion scheme can be generalized to give mode- 
mode coupling equations with the bare interactions replaced by matrix elements of 
the Liouville operator. 

In our illustrations we have been concerned with spin systems. But the given for- 
malism also applies to Fermi or Bose systems, if for the statics of Fermion systems 
the fluctuation-dissipation theorems are used in an adequate way. Thus, the possibil- 
ity arises to relate expansions in Mori's theory to other many-body treatments, and 
especially find a link between projection-operator expansions and those in Green's 
function theories. 

Summarizing we think that the main advantage of the presented expansions is, that 
they lead to approximations for the Heisenberg dynamics. These can be used to derive 
equations of motion for expectation values which go beyond linear response. 



A Second order approximations to the residual forces 

A.l Expansion into powers of e ~ with fixed correlations E(t) 

The expansion of the residual forces in powers of e, with fixed dynamic and static 
correlations, is obtained from (2.17). The first order result was given in (2.20). Up 
to second order one finds 
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f(t)=tQV^{GE,GE} 

+ iQV {i) {iQV {2) {GE, GE} <g> S, GE} + *Q\/ (2) {G~, *Q\/ (2) {G~, GE} <g> 2} 
-2g\/ (2) {GS,G'S}®(G'|G')- 1 (G'K\/ (2) {G'S,G'S}) + 0( e 3 ) . (A.l) 

The approximate /(t) now is an element of a linear space spanned by the products 
QGG and QGGG. The selfconsistent equations for the dynamic correlations E(t) are 
given by Eq.(2.9) with the force correlation functions (f\f(t)) calculated from (A.l). 
A general feature of the expansion (A.l) of f(t) is that terms of second, and higher 
order, will contribute to the projections onto the forces at t — 

f = iQV^{G,G} . (A.2) 

That means they will contribute to the force correlation functions. In the sense of an 
orthogonal decomposition of /(£), the coefficients of QGG are changed by each order 
of the iteration. This is, why it seems to be reasonable, to sum up all contributions of 
higher order terms which are parallel to the products QGG. To make this reasoning 
explicit, we introduce a basis F v for the QGG by 

Fu = QJ2 oc v ,AG,G x - (G,G X )) (A.3) 

which is abbreviated by 

F = Qa{G, G} - (a{G, G}) . (A.4) 

Then f(t) can be decomposed as 

f(t)=F(F\F)-\F\f(t))+Qf(t) , (A.5) 

where 

Q = {(1 - \F)(F\F)-\F\}Q , (A.6) 

and (F\f(t)) is no longer affected by the approximations of Qf(t) which are orthogo- 
nal to (A.3). In the next section A.2 we will calculate f(t) to second order with fixed 
E(t) and fixed components (F\f(t)). 



32 



A. 2 Expansion into powers of e ~ with fixed correlations H(t) and T(t) 

For a treatment of f(t) according to (A. 5) it is expedient to decompose fit) into a 
sum 

f{t) = e iQLQt iQLG = K(t)iQ^ , (A.7) 

where 

K(t) = e iQLQt F , (A.8) 

and the matrix is defined by 

fi (r) = (F\F)-\F\LG) = (F\F)-\F\V {2) {G,G}) . (A.9) 

Here use had been made of the fact thatf^] 

QLG = Q(V {1) G + V (2) {G,G}) =FVl {r) . (A.10) 

Then corresponding to (A. 5) we can write 

K(t)=FT(t) + QK(t) (A.ll) 

with 

T(t) = (F\F)- l (F\K(t)) = (F\Fy\F\e lQLQt F) . (A.12) 

In the sequel we will focus on the expansion of the quantity K{t), which is of course 
equivalent to the expansion of the residual force f{t). 

For K(t) we will set up an exact equation with fixed H(t) and T(t) which can be 
iterated. As a first step we cast K(t) into the form 

K(t) = FT(t) + Qe lLt QiLF ® T(t) - iQe iLt F ® (F\F)-\LF\K(t)) (A.13) 

which can be verified by Laplace transform and the completeness relation 

F(F\F)-\F\QGG) = QGG - (GG) . (A.14) 
14 From (l\LG) = and (l\G) = it follows (1\V^{G, G}) = 0. 
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Next we take the Heisenberg dynamics of G and F, which can be expressed as 

G(t)=GE(t) + K(t)iQ (r) ®H(t) (A.15) 
F(t)=G~(t) ®itt {l) T(t) + K(t) (l - ®n (r) Z(t)rt l) <g>r(*)) . (A.16) 

Here denotes the matrix 

SI® = (G\G)- 1 {G\LF) = (G\G)- 1 (VW{G,G}\F) , (A.17) 

and use has been made of Eqs.(2.9) and (A. 7), and 

G = F{t)itt {r) +G{t)iQ = GE{t) + f{t) + f{t)®E{t) . (A.18) 

Calculating LF from the definition (A. 3) and inserting Eqs.(A.15) and (A.16) into 
the right hand side of Eq.(A.13) one arrives at an exact system of nonlinear equations 
for K(t) 

K(t) = FT(t) (A.19) 
+ Qa s {-GE(t) <g> Q^T{t)Q^ + K{t)iQ^ (l - <g>5(*) <g> Q {l) T{t)Q^) , 

GE(t) + K(t)iQ {r) <g> E(t)} <g> r(t) 
- QK{t)iQ {r) <g> S(t)(G|G) _1 (GMFin (r) , G}) <g> 
+ Qtf(i) (l - ®ft (r) H(t)fi w <g> r(t)) ® (Fli 71 ) -1 (a s {™ (r) , 

where the abbreviation 

a s {AB} := a{A 5J B} + a{ J B,A} (A.20) 

has been introduced. This system corresponds to Eq.(2.17) for f(t). The difference is 
that in Eq.(A.19) the correlation functions E(t) of the observables and the correlations 
T(t) (A. 12) of the residual forces (A. 8) appear. 

The system (A.19) now can be iterated. Remembering that according to (A. 9) and 
(A.17) fi (r) and fiW are of first order in V®, we find from (A.19) for fixed E(t), T(t) 
(in treating T(t) fixed the summation of higher order contributions is performed): 

K{t) = FT{t) + Qa s {FT{t)iQ (r) ,GE{t)} ®F(t) + C(e 2 ) . (A.21) 

According to Eq.(A.7) it is a second order result for f(t). The main point is, that 
the approximation (A.21) conserves the exact relation 
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(G\G)~ 1 (f\f(t)) = (G\G)-\iQLG\F)(F\F)-\F\K(t))in^ 

= fi ( °r(t)fi (r) . (A.22) 

The matrix T(t) can be calculated from the differential equation 

f (f) = i(F\F)-\F\LF)r(t) - {F\F)-\iLF\QK{t)) (A.23) 

which follows from Eqs.(A.ll) and (A. 12). Inserting the approximation (A. 21) for 
QK(t) and using the definition (A. 4) we find 

t(t)=i(F\F)-\F\LF)r(t) 

-(F\F)- 1 (a s {Firt r \G}\Q(a s {FT(t)irt r \GE(t)})) ®T(t) . (A.24) 

Together with the dynamic equation for E(t) 

E(t) = E(t)iQ - E(t) <g> fi«r(t)fi (r) (A.25) 

we have a coupled system for the correlations F(t) and E(t). As Eq.(A.25) can be 
solved by Laplace transform in terms of T(s), Eq.(A.24) can be viewed as a nonlinear 
integral equation for the Laplace transforms T(s). It should be noted that the ap- 
proximations (A.24), (A.25) have the exact second order moment of T(t) or the exact 
forth moment for E(t) respectively. For later use we note the Heisenberg dynamics 
for G(t) which result from the approximation (A. 21). One obtains 

G(t) = GE(t) + f(t) ® E(t) = GE(t) + FF(t)iQ^ ® E(t) 

+ Qa s {FT{t)iQ {r) ,GE(t)} ®T(t)^ {r) ®E(t) . (A.26) 

A. 3 Expansion of the residual forces f(t) 

We want to show that the results of section A. 2 can be obtained in a simple way by 
extending the set of observables G to G and applying the procedure of section 2.2 to 
/. Let us choose a set G by 

G = {G l/ ,F ll } (A.27) 

where F u is defined by (A. 3). Then we have a projection operator (A. 6) 

Q = 1-P=1-|G)(G|G')- 1 (G'|-|F)(F|F)- 1 (F| . (A.28) 
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The matrix can be written in block matrices as 



-'GG ^GF 
-FG ^FF 



(A.29) 



Similar the frequency matrix Q reads 



n 



^GG &<GF 

Qfg &<ff 



(A.30) 



where £l GG — and CIgf and Qfg coincide with the definitions (A. 17) and (A. 9) 
for fiW and fiM. 

First we study the derivatives LG and prove that the conditions (2.14) of section 2.2 
apply. The derivative 



LG = V {1) G + V (2) {G,G} 



(AM) 



is a linear combination of G,, and F„. Hence in abbreviated notation we have 



LG — GVl GG + FVtpQ 



(A.32) 



The derivatives LF are found from (A. 4). Together with (A.32) one yields^] 



LF = G [(G\G)- l (G\a s {Gn GG , G}) - ^(GI^-^G^jG, G}) 

+ F [(F\F)-\F\a s {Gn GG , G}) - Q GG (G\G)-\G\a{G, G}) 

+ a s {Fn FG ,G} . (A.33) 

Eqs.(A.32) and (A.33) together show that LG has the desired structure 

LG = V (1) G + V^{G,G} , (A.34) 



where the matrix elements of V^ 1 -* and can be taken from (A.32) and (A.33). A 
special feature of is, that it just couples F and G, but no couplings FF occur. 



15 We take the case that for symmetry reasons the constant term in (A.33) vanishes. If it 
does not vanish it can be handled without difficulty. 
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Having shown that the condition (1.1) holds, we use Eq.(2.17) to expand the residual 
forces /. As 

V (2) ~ Q FG ~ V (2) ~ e , (A.35) 

the interaction still is of first order in e. Expanding / we want to keep fixed 
E(t) and F(t) in the sense of section A. 2. Therefore we must express the correlations 
E(t) in terms of E(t) and T(t). From the definitions of E and (A. 15) and (A. 16) one 
obtains 

E GG (t)=E(t) 

E FG (t) = r(t)inw ®E(t) 
z GF (t) = z(t)®inWr(t) 

E FF (t)=T(t) (l - ®n (r) Z(t)tt (l) <g>r(t)) . (A.36) 

This means that fixing S(t) and T(t) is equivalent to fixing the matrices E GG , E GF , E FG , 
but due to Q( t \Q^ in (A.36), the submatrices E FG and E GF are to be treated as 
first order terms, whereas E FF has a zeroth and a second order contribution. Hence 
expanding / in Eq.(2.21) up to first order yields 

/ = QV^{(GE)°, (GE) } + 0(e 2 ) (A.37) 



where 

(GE) G = GE GG , (GE) F = FT(t) . (A.38) 
Explicitly the residual forces (A.37) read 

f= (/f) = {Qa s {FT(t)iQ FG ,GE GG (t)}) ' (A ' 39) 
The selfconsistent equations for E(t) are found from and (f\f(t)) to be 



-gg = E GG i£l GG + E GF iVt FG 

*->gf = E GG iQ GF + E GF iVt FF — E GF <g> 7i?i? 

^fg = E FG iVt GG + E FF iVt FG 

Iff = E FG iQ GF + E FF iQ FF — E FF <g> 7^ (A. 40) 



37 



where 



lFF {t) = (FIF)- 1 (iLF\Qa s {Fr(t)in FG ,GZ GG (t)}) ®T(t) 



(A.41) 



The approximation for the Heisenberg dynamics for G(t) are obtained from 



G(t) = GE GG (t) + FE FG (t) + Qa s {FT(t)tQ FG , GE GG (t)} ® E FG (t) . (A.43) 



Now the results of sections A. 2 and A. 3 can be compared. One sees that the Heisen- 
berg dynamics (A. 26) of section A. 2 agree with (A.43) of section A. 3: Inserting the 
relations (A. 36) into (A. 26) one obtains the results (A.43). Furthermore, using the 
relations (A. 36) as a definition for the left hand sides, one can convert the equations 
of motion (A. 24) and (A. 25) into the Eqs.(A.40) for E. Thus we have demonstrated 
that the approaches of sections A. 2 and A. 3 are equivalent. It is clear, however, that 
one can also expand / with E(t) being fixed and taking E GF , E FG of first order, which 
replaces T(t) in (A. 38) by E FF . Alternatively, one can directly use the expansion of 
section 2.2 leading to / = QV^ 2 \GE,GE), c.f. Eq.(A.37). Such expansions sum up 
further higher order contributions and would include the results of section A. 2. The 
spin model of section 3.1 is treated in this way. 

We conclude this section making some remarks upon the static correlations which 
occur in the treatment with the extended set G. Let the approximate / be a lin- 
ear combination of Q(GF) (e.g. Eq.(A.39)), and take the commutator forms for Q 
and (f\f(t)), then the correlation functions E(t) depend on (G\G) and (GG), and 
(G\GG), similar to the case of Eq.(2.35). The essential point however is, that ow- 
ing to the commutators [G, G], (G\GG) only enter as combinations (G\GF) and 
(F\GF). Treating the correlations along the lines of section 2., with exact equations 
relating (G\G), (GG) to (G\GG), (G\GF), (F\GF), and using the approximations for 
(G\GG), (G\GF), (F\GF), one sees that the system is closed, if (FGF) is included. 



B n— th order approximation and exact solution for the spin model of 
section 3.1 



G(t) = GE(t) + f(t)®E(t) 



(A.42) 



to yield 
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B. 1 Self consistent equations for 



Take the orthogonal set G = {G , G q ,i, . . . , G q , n ', 9^0} defined by (3.31), then the 

— * 

time derivatives of G q>u read 



iLG = 

iLGg^l = G q ,2 

iLG qiV = -c v G q ,v-i + G q ,v+i v = 2, . . . , n - 1 

— * — » 3 ( — *—»—»— » \ 

iLG q ^ n = —c n G qtn -i + — ^= (^Go x G 9jTl — G 9jTl x Go J (B-l) 

with c v defined by Eq.(3.39). This means that in the interaction 

LG = V (1) G + V (2) {G,G} (B.2) 

just matrix elements V^fin an d Ki,no are different from zero. Therefore the approxi- 
mation (2.20) with e = J for the residual forces yields 



/,» = !/ = l,...,n-l 

fl n {t) = Q^(G x G,, B - G q , n x G ) Q S nn (i) + 0(e 2 ) 

= Gl n+1 E nn (t) + 0(e 2 ), a = xoiy,z (B.3) 

leading to the memory matrix 

(f a |p (t)) 

-MO = XX = ^n^c„ +1 S nn (t) + 0(e 3 ) . (B.4) 



The equations of motion of the matrix E Ufl (t) (2.9) follow from (2.7) with 

iQ. vll = -CpS^-x + S^+i (B.5) 



and (B.4) for 7^(t). In the approximation for the Heisenberg dynamics (2.23) the 
matrix elements E u i, v = 1, . . . ,n and E nn enter. We therefore regard the Laplace 
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transforms of S„i and S n „ which can be written as 



sS 2 i + C3S31 



—11 



= 



(B.6) 



sE vl + c u+1 E u+hl - E v - ly i = v = 2, . . . , n - 1 



and 



= 

+ c 2 H„i = 



^2 , _ 1 



ran— 1 



(B.7) 



i/ = 2,...,n-l 



Eqs.(B.6) and (B.7) provide the selfconsistent solutions. The coefficients Cj, can be 
found from the orthogonality of the G qi „ to yield 



ci = l 



2 ' 

¥5-1 



^0 = 1, ¥>i/ 



1/ > 1 

CKi «2 ••• Ct u 
Q>2 «3 • • • Ot v+ i 



OL2v = 

"2^+1 = (-l) u m 2u 
v > 0, 



«j, C^+i . . . OLi v -\ 

where mi v denotes the exact 2u-th moment 

"fa I / T \2v r~iot \ 



m 2v = 



(B.8) 



(B.9) 



which can be expressed in the thermodynamic limit N — > 00 by the second moment 
(see (B.29)) 

m2 ^K)^^ • (B.10) 
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B.2 Explicit Solution o/E< n )(s) 



The system (B.6) can be solved in terms of S n = E^ n \ With help of the polynominals 
A u (s) and B u (s) defined by Eq.(3.38) one directly finds the result (3.35). For the 
solution of the system (B.7) we first take the equations for v — 1, . . . , n — 1 expressing 
^nu+i by S„i which gives 

E nu+1 = B u E nl u = l,...,n-l . (B.ll) 



With help of (B.ll) and (3.38) and 



'nn— 1 



B 
B 



n-2. 



(B.12) 



n-l 



the last equation of (B.7) can be converted into a quadratic equation for E nn 

Cn+li^nn) 2 + D ~ ^nn -1 = , (B.13) 
-Dn-1 



which has the solution 



B ri 



B n -i \ 



B r , 



B 



+ 4c, 



n+l 



n-l. 



2c n 



+1 



(B.14) 



The other solution is not a Laplace transform of a function E nn (t) regular at t — 0. 
Combining Eq.(3.35) ior v = n with (B.ll) for z/ = n — 1 gives the solution (3.37) of 
section 3.1 for S^(s) which can also be written in terms of E nn as 



- D n ~r t -n+l^nn £l n-l 

if, following from Eq.(3.38), the relation 



(B.15) 



A n B n -i - A n _ x B n = (-l) n+1 Cl ■ ■ ■ c n (B.16) 



is used. 

The analytic properties of S(")(s) follow from S nn (s) and the representation (B.15). 
First one proves that E nn (s) is holomorphic for Res > 0, and it holds 

ReZ nn (s) > for Res > . (B.17) 
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To see this, we take the definitions of the polynominals B v (s) writing 



^ = S + ^ r , c„>0 (B.18) 

D v— 1 -5 



which by induction leads to the result 

B 



Re— — > 0, Res > . (B.19) 
B v -\ 



Together with the property, that the B v (s) have no zeros for Res > [9], the mapping 
(B.14) of B n /B n _i onto E nn completes the proof. For S^(s) in Eq.(B.15) we use the 
recurrence relations (3.38) writing 



B n + c n+ i^ nn B n -i B n -i + s+c ^" 1 ^ nn B n ^2 



(B.20) 



where from (B.17) it follows that 



Re 1 - > for Res > . (B.21) 

S + ^.(-Il^^^SJ 



Iterating this procedure one sees that (B.17) induces E^^s) to be holomorphic for 
Res > and to have the property 

ReS (n) (s)>0 for Res > . (B.22) 



As a last point it is not difficult to see that for s approaching the imaginary axis, 
Res — > +0, the function S(s) keeps to be finite. 



B.3 Exact solution E(t) 



To derive the exact solution for the dynamic correlation function 

TrG a „G°(t) 



42 



we start with an exact differential equation of second order for the Heisenberg dy- 
namics of G q (t) 



In the Hilbertspace of our model (N spins 1/2) we use now a partition of the iden- 
tity operator into operators Ps projecting onto the eigenspaces of So ■ So leading 
immediately to 



' J " 2 



+ 4l-=l S(S + l)TrG a _ q P s G a q (t) = (B.25) 



with the solution 



TrG°L q P s G°(t) 

+ (^T^V^ + ^y^G;(0))e- 2 «^ . (B.26) 
By calculating the traces which are independent of q and a = x or y, z, with 

^-I^iLL 1 ■ (R27) 



and evaluating the sums we arrive at 



N-2 



E im ( t ) = NU-Lt) (E±l^('t)-i) (B.28) 



or, taking the thermodynamic limit N — > oo 



E(t) = lim S (JV) (f) = e- j2 * 2 / 2 (l - JH 2 ) . (B.29) 

A<— >oo 
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C Dynamic and static correlations of a Heisenberg ferromagnet 



C.l Self consistent equations in lowest order 



By the approximation discussed in section 3.2 the dynamic correlation functions 
S*(t) and S ± (t) of a Heisenberg ferromagnet depend on the following static correla- 
tions (5S Z \5S Z ), (6S Z 5S Z ), (5*15*), (5"5 Z ), (5S Z \S-S+), (SWISS'S*), (S-S + \5S Z ), 
(5*15-5+^) (see Eqs. (3.69)-(3.78)). 

In this appendix we want to point out now, that one can close the dynamic equations 
(3.69)-(3.75) of section 3.2 by a set of relations determining the static correlations 
in terms of the dynamic ones. The procedure is analogous to the the general scheme 
discussed in section 2.5, the only difference is, that the iteration with Vq +eV x has 
led to residual forces, which have contributions of the form GGG. Therefore there 
appear higher-product static correlations, e.g. (S ± \S~S + S ± ), and the set of static 
quantities has to be extended in order to find a closed system of equations. 

Again we start with exact relations between (G\G) and (G\GG) on the one hand, 
and between (G\G) and {GrG) on the other hand (compare Eqs. (2.36) and (2.37) of 
section 2.5)0 



(ss!\ssf) = -^^^2(5Sf\S2SU) - j2Tw^ {ss ' lss ' 2SS '- 2) (C ' 2) 



N = 2P(J - J^SflSf) - e^(J!_ 2 - J 2 )(St\SS z St 2 ) (C.3) 

iV 2 

= (5S Z \5S Z ) J du-J^-^liu) (C.4) 

(51,5+) = (5 + |5 + ) [duJ^-Et(u) . (C.5) 
J eP" — 1 



16 For the longitudinal case the relation analogous to Eq. (C.3) would read 

= £(Ji- 2 - J 2 )(SS z 1 \S 2 Si_ 2 ) (C.l) 

2 

and would not supply any further information, because by symmetry Eq.(C.l) is identically 
fulfilled. We therefore have replaced Eq.(C.l) by (C.2), following from 

^(<5^|5 2 -5 1 _ 2 ) = 
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Using now the approximation (3.62) and (3.63) in Eqs.(3.65) and (3.66) of section 
3.2 for the Heisenberg dynamics of 5S z (u) and S + (u) 



2? - - 

ss^u) = ss^i(u) - Tr Y,Q s 2SU(Ji-2 - J 2 )(0i2 ® SDH (ce) 



S+H = S+S+( w ) + £( J x _ 2 - J 2 )Q {<5SJS+ 2 ((S^S+_ 2 ) ® S+) H- 

- |E^(4"4 + - 3 )^ + - 2 (J 2 -3 - J 3 ) [((023 ® S*) S+_ 2 ) ® S+] (c)| (Q.7) 
the equations analogous to Eqs.(2.28) and (2.29) of section 2.5 explicitly read 

(A\5St)=a 1 ([S z 1 ,A^]) + j:b 12 ([S 2 St 2 ,A^ (C.8) 

2 

(A|S+)= Cl <[S+ 

+ e^d 12 ([^+ 2 ,At])+6^e 123 ([4"^i + -2^]) (C9) 

2 2,3 

(^)=a 1 ([^,A]) + 5:6 12 ([4-5 1 t 2 ,A]) (CIO) 

2 

(A^>=c 1 ([5?,A]> 

+ e^d; 2 ([^+ 2 ,A]) + e^e 123 ([4-4 + - 3 ^ + _ 2 ,A]) (C.ll) 

2 2,3 

with coefficients 



al= y^ 1H -? (*7fl*7f) 612 (c - 12) 

612 = -|(^l-2 " J 2 ) / ^(012 ® HDH (C13) 

« > -%^Wf d » + % - twr 2> e - (ai4) 

ei 23 = - (f)Vi-i - Ji)(J 2 - 3 - J 3 ) / ^ ® S83+_ 2 ) ® H^} ( W )(C,16) 

and analogous expressions for the coefficients 01, b 12 , Ci, d 12 and e 123 . The only differ- 
ence is, that the denominators f3u in Eqs.(C.12)-(C16) are replaced by (e^ w — 1). 
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We now usef 17 ] Eq. (C.8) for 



A — Si_ 2 ^2 



(C.17) 



Eq. (C.9) for 



A 



Q+ Q— cH 
°l-2°2-3°3 



(C.l* 



Eq. (CIO) for 



A 



^-3^lt-2 < ^'2-l 
, >J-4' : '4-3*- , 3-2'- , 2-l 



(C.19) 



and last not least Eq. (C.ll) for 



(C.20) 



and evaluate the commutators on the right hand side with help of 



S\ > $2 



NS^ 2 + eSSt 



3=0 1+2 



(C.21) 
(C.22) 



Then it is easy to see, that Eqs.(C.2)-(C.5) together with Eqs.(C.17)-(C.20) con- 
stitute a system of 12 equations determing the 12 static correlations (5S Z \SS Z ), 
(5S Z 5S Z ), (S-\S-), (S-S+), (5S Z \S-S + ), (S-\5S Z S-), (S-S + 5S Z ), (S-\S-S+S-), 
(5S Z 5S Z SS Z ), {S-S+S-S+), (S + S + SS Z 5S Z ), (S~ S+S' S + 5S Z ) in terms of the dy- 
namic correlation functions H 2 (t) and H ± (t). Thus together with the equations of 
motion (see Eqs.(3.69) and (3.70) of section 3.2) 



-7i ® si 



(C.23) 
(C.24) 



17 For the transverse case it suffices to regard correlations like (S 1 \ . . . S 1 ). The correlations 
(S>i\ . . . Sf ) follow from symmetry arguments, e.g. (Sf \Sf) = (5^~|S^~). 
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with 7I and 7^ given by Eqs.(3.74) and (3.75) of section 3.2, we have found the 
desired closed system of equations coupling static and dynamic correlations. 



C.2 Zeroth order solution 



The easiest way to find the zeroth order contributions in e to the dynamic- and 
static correlations of a Heisenberg ferromagnet is to use perturbation theory to the 
Heisenberg equation of motion (3.58) and (3.59) of section 3.2 for 8S z (t) and S^t). 
For demonstration purposes we want to sketch, how these simplest approximations to 
the correlations can be found from the zeroth order solution of the rather complicated 
coupled system of equations constructed in appendix C.l 

First the zeroth order of the frequency uof^ (see Eq.(3.73) of section 3.2) can be found 
from Eq.(C.3) to yield 

^! 0) = - , ; = 2(J - Ji) =■ &! . (C.25) 



Since the correlation function 7~ L (r) of the transverse residual force is of second order 
7 _l (t) ~ e 2 , Eq.(3.75) can easily be integrated 

~±(°)( r ) = e w (C.26) 



by Eq.(C.5) immediately leading to 

^(^ 1 + )(°) = (e' 5 --l) _1 =:n 1 . (C.27) 



To find the zeroth order of the longitudinal correlations, we first notice, that (C.2) 
leads to 

(SSilSSt)® = ~*E(6Sl\S;St 2 ) l0) (C.28) 



and (C.17) to 

(SS^SU) (0) = N (aS 0) - 6ff_ 2 ) (n 2 - ni_ 2 ) (C.29) 

with 
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fll -^y^" 1 (55f|(55f)(o) 612 (C - 30) 

&S? = (a) 2 " ^1-2) Pr/ i (0 12 ® Sf)^ ( W ) . (C.31) 

To evaluate the principal value integrals they are expressed by 

p I [fo*( u ) = *(* > (r l )-* < (-T l )), rj^+o (c.32) 

J CO A 

with the abbreviation 

00 

*>/<(s) := J dTe TST ^(±r) . (C.33) 



The expressions S>?<(s) and (0i2®S^)^ < (s) can be found from the defining Eq.(3.69) 
together with (3.74) of section 3.2 using the zeroth order result (C.26) of Sf (r). They 
depend on the static correlations (6S*\S 2 S?_ 2 )^ and ((SS^S*)* ). Thus Eqs.(C28) 
and (C.29) together with the definitions (C.30) and (C.31) lead to a coupled system 
of equations, which can be solved by using Eq.(C.32) to yield 



(5S*\S 2 SUr=N ^ n !~\ (C.34) 

(5S Z 1 \5S Z 1 )^ = -T i n2 ~ ni ~ 2 . (C.35) 
2 P(u 2 - W1-2) 

This implies 

S?° >( T ) = f E ^ 2 ~ ni ' 2 1 1 E ~ n2 ~ ni " 2 e'^-^-^ . (C.36) 
V2 (3(^2 -0J1-2)) 2 /3(w 2 -wi_ 2 ) 



The results for the expectation values 

(<SSV<S?) (0) =5> 2 (l + ni_ 2 ) (C.37) 
2 

(SZ 2 StiSSZ)W = -Nn 2 (l + m_ 2 ) (C.38) 
follow from Eqs.(C.4) and (C.19) with Eq.(C36). 
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